NACA TN No. 1443 


J . /Y iu /o' ^ 


GOVT. DOC. 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 

No. 1443 

SHEAR LAG IN A PLYWOOD SHEET -STRINGER COMBINATION 
USED FOR THE CHORD MEMBER OF A BOX BEAM 
By Palamede Borsari and Ai- ting Yu 
Massachusetts Institute of Technology 


Washington 
March 1948 


MAR 31 1948 


'business, suencl 

& technology DEPT 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO. lU+3 

SHEAR LAG IN A PLYWOOD SHEET -STRINGER COMBINATION 
USED FOR THE CHORD MEMBER OF A BOX BEAM 
By Palamede Borsari and Ai-ting Yu 

SUMMARY 


Theoretical and experimental investigations were made of the 
distribution of strains in a plywood sheet-stringer combination used 
as the chord member of a box beam acted upon by bending loads. The 
theoretical solution was obtained with the help of the principle of 
minimum potential energy and certain simplifying assumptions. Strain 
measurements were made on a built-up box beam by means of electrical- 
resistance strain gages connected with strain indicators. A very 
satisfactory agreement between the theoretical and experimental strains 
was obtained. 


INTRODUCTION 


The development of monocoque and eemimonocoque structures in 
airplanes has introduced the shear-lag problem, which, in brief, results 
from the fact that box beams with thin, wide chord members do not follow 
elementary beam theory. The assumption of uniform normal stress distri- 
bution on fibers equidistant from the neutral axis of any transverse 
section is not justified for these beams* The deviation from uniform 
normal stress is caused by the shear deformation in the chord member 
of this type of structure, an effect not provided for in simple beam 
theory. 

Various theoretical treatments have been carried out and numerous 
experiments have been conducted to check the adequacy of them, especially 
on isotropic or orthotropic metal structures. Very little information, 
however, is available on plywood construction. For plywood sheet -stringer 
combinations there appear to be no data. 

The present investigation is a condensed version of a thesis 
presented to the Massachusetts Institute of Technology in partial fulfill- 
ment of the requirements for the degree of Master of Science in Aero- 
nautical Engineering. The theoretical development presented in this 
thesis is due to Dr. Eric Reissner, who outlined to the authors the 
several steps necessary to obtain the solution. 
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SYMBOLS 


x 


y 

z 


L 

2d 

h 

t 

b 


•^Bt 

E X 

E st 

E.. 




u sheet> 


^eb 

w web 

M x 



€ web 


xy 


spanwise rectangular coordinate measured from tip to root in 
midplane of cover sheet 

transverse rectangular coordinate measured from longitudinal 
plane of symmetry to edge in midplane of cover sheet 

coordinate axis perpendicular to neutral plane of beam measured 
from bottom to top cover 

beam span 

beam width 

beam depth 

cover-sheet thickness 
side-web width 

area of cross section of nth stringer 

modulus of .elasticity of cover sheet in x-direction 

modulus of elasticity of nth stringer 

modulus of elasticity of side web 

modulus of rigidity of cover sheet referred to shearing 
stress acting parallel to x- and y-axes 

v displacements in midplane of cover sheet in x- and y-directions, 
respectively 

displacement of side web in x-direction 
displacement of side web in z-direction 
applied bending moment 
strain components in cover sheet 
longitudinal strain in stringers 
longitudinal strain in side webs 
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DEVELOPMENT OF THEORY 


By following the suggestions given by Reissner (reference l), 
a complete solution of the shear-lag problem for the case of a sheet- 
stringer combination used as a chord member of a singly symmetrical box 
beam of constant cross section is presented. Further generalization of 
the method may be carried out with little difficulty, since the general 
procedure remains essentially the same for any caramon shape of beam 
or load. 


Principle of Method 

The method consists in expressing the total energy of the system 
as the sum of the energy of deformation (the volume integral of the 
strain-energy function) and the energy of external forces- Thus, 


n = + n e 


( 1 ) 


Following this, the differential equations of equilibrium are obtained 
from the theorem of minimum energy stated as: The displacement which 

satisfies both the differential equations of equilibrium, as well as the 
conditions at the bounding surface, yields a smaller value for the potential 
energy of deformation than any other displacement, which satisfies the same 
conditions at the bounding surface'' (reference 2) . This can be written 
as 


Sit = 0 


( 2 ) 


Assumptions for Simplification of Problem 

The expressions for and n e are based on the following simplifying 

assumptions: 

(1) It is customary, in dealing with problems of this nature, to 

neglect the stretching of the sheet in the transverse direction, 

since it is very small as compared with the stretching in the spanwise 
direction. This assumption has been used even in "exact" analytical 
solutions, as shown in reference 3* 

(2) Since if = 0, v would be a function of x only. 

Here, however, we assume that the displacement 


v = 0 


( 3 ) 
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(3) Instead of assuming a parabolic distribution of stress as in 
reference 4, Reissner suggested a parabolic distribution of displacement 
of the sheet in the transverse direction. These displacements can be 
written in the form 


(4) The linear side-veb displacements are determined in such a way 
that cover-sheet and side-web normal displacements coincide along the 
line of junction, that is. 


(5) The sheet iB free from shear parallel to the z-axis. Such shear 
is carried by the webs only. 

Equation (2) furnishes the differential equations for the determination 
of the unknown functions uq, U]_, and U2 of equations (4) and (5). 

Since 



(«0 



(5) 




Su dv 
= bj + c>x 


assumptions (l) and (2) are expressed by 




and 
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Hence, 



— 


( 6 ) 




J 


Then, since Uq, u^, and. ug are functions of x only, equations (3) 
and (4) imply a parabolic distribution of normal strains 


A cantilever box beam (or an equivalent simple beam having double 
the length of the cantilever span) is considered. The cross section 
is constant along the span and the beam is acted on by a given distri- 
bution of bending moment. The box beam is singly symmetrical. It has 
two side webs, with a sheet -stringer cover sheet on one side. (See 
figs. 1 and 2.) The centroids of n.1 1 the stringers are assumed to be on the 
midline of the sheet; hence bending effect on the stringers is neglected. 

The potential energy of the system is obtained from 



(7) 


and a linear distribution of shearing strain 



(6) 


along any transverse section. 


Formulation of Problem 



( 9 ) 
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^web 


P L 

Jo 


h 

bK^ e v 2 dz dx 

Jo 


( 10 ) 


K 

*st = 5Z 
n=l 


n L 


A st E st e st d* 


( 11 ) 


Neglecting the shear in the side webs yields 


7xz web 


Su veb ^eb 

+ 

dz dx 


^seb = jKeb 

dx dz 


From equation ( 5) , 


^ u web 

dz 



Therefore, 


dv 


veb - -&*> 


dx 


Again, 


_ f M d^v 

rt load I , 2 ^ 

Jo dx* 
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Therefore, 


load 


FL — ' l ~ dx 
^ dx h 


( 12 ) 


Finally, 


1 

1 L 

d 

(Z 2 

1 L 

n - 2 

U 

0 i 

t 

-d 

(Ex € x + dx dy + 

o Jc 


hE w € w dz dx 


* Z f A st\t‘st 2 to - f “xSTTito 


n=l JO 


(13) 


For minimum energy condition, fham equation (2) , 
L n d 


bn = 


l J t ^x de x + ^xy 8 ^xy) & 


L fl h 

+ I | 2hE w e(x,w)6e(x,w) dz dx 

'0 *0 


Z f 2A 8t E 8t € St 5e st 
n=l a 0 




( 14 ) 
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From equations (6), (7), and (8), the variation of strains is given by 


6e x 

sheet 



+ 




6 



5 ^xy 



&e 


■Sreb 





Theref ore. 



8uq 

sheet dx 



+ 


f y^i Su O . /^ &u 2\ 








Substituting these relations into equation (14) yields 



_Bx Vd^y d 2 y dx Vdxy 



L „ h 


2b E* 


ilo do 


s 2 <g) <1 - J) af 3 CgKs ' i) 




a • *>' s> .©; 


dz dx 


K 

z, 

n=l u 


2 A 0 4. E 


st -^st 


^0 d /^o\ . 

dx \3xy 



d /^£) + (\ . y jL 
dx \dxs \L d^y dx \3xy 


y n 2 \ dug /SuA 
+ " d^y dx * vdx ) 


dx 


f L 

Jo 


*x 


/ 6 Ut\ , 

d ( - — ) ■* dx = 0 

\dxy h 
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Integration with respect to y in the first integral and with respect 
to z in the second integral can now he carried out. 

1 d 

dy = 2d 

J -d 
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Substituting these values and grouping the terms accordingly, yields 




2 

3 


n 

t Ej d + 

1 







4t 

+ — Gup 
3<1 


6up 


b Evh du 0 . /^ u l\ 

“2~ sr b \&y 


b Eyh /^ 6u 0\ 

~ ST d v^7 
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Let 


n 

<P| = t E x d -t- ~b Eyh + y~ A st E st 

1 



Integrating by parte gives 


c)uq Sug dvu 

^ aT 6u ° + 92 aT 8u2 + *3 aT 5u o + ^ aT 8u 2 


b E^h 


/^Uq ^ 2 ^ u i ~\ , 

sr 5u i • sr 8u o + 3 sr 5u y + 4 ^ 5u i 


-I 


J 


a 2 - 


a 2 - 


a 2 * 


a 2 i 


iiq ° Uo 0~i*2 o“u2 

*Pl TIT 5u O + <fe 7 p~ 5u 2 + «h “o" 8 ^0 + 8u 2 

ox j ax* ax* 


4t _ 

3d GU2 5u 2 


b Eyh < 


^ 2 up a 2 ^ a 2 ui 'N 1 
6u i ' p- 5u o + f S'i) -as 8u i 


dx = 0 


The first bracket is equal to zero from the following assigned boundary conditions: At x = L, 


and at x = 0, 


Suq = &u-j_ = 6ug = 0 


a Ul aug 

ax ax ax ^ 0 
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The integral thus becomes 


1 L 

Jo 


f ^ 2 uq d 2 u 2 b Evh a 2 uj\ 

+ % ' 2 *£) 5u ° 


2 2 
h EJi d uo t Eyh d u x ! 

+ 2 dx 2 + 3 ax 2 2h ZiJ 1 


*2 


^Uq 


S 2 


dx 2 + ^ ax 2 3d 


4t 


rr Gu 2 ) &U2 


dx = 0 


Since bu^, 6u2, and &uq are arbitrary. 


' 2 % b Eyb ° 2u l 

"*• * *3 a? 2 SF ‘ 0 


A> d ‘ 


2 2 

* Eyh ^ U Q + b Eyh a u x _ i_ = 0 


dx‘ 


3 ax 2 " 2h 5x 


^2 
o u r 


"2 4t 


"fe S? - + ^ £T ' 3d GU2 = 0 


Let 


2 <?L 
b Eyh 


X 2 = 


2q>} 


b l^h 
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Then 


By eliminating 


= 


3 dqfe 

4t G 


H = 


3 dtp4 
4t G 


4 * = -■ 


* E v h 


X 1 




bx 


b 2 u 2 b \ 

+ X2 ^-^--° 


2 d U 1 dM 0 \ 

■— + , + x 5 s = 0 ? 


3 


X3 + 57 ” U2 = 0 




uq, the following equation is obtained: 


(15) 



(16) 
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Again, "by eliminating Uj_, the following equation is obtained.: 




% ■ *) I 


^2^ ^ 


O Ug 

ax 2 


3 l l • i) ”2 ■ X 5 X 3 s 


By introducing the parameters and 0g 83 follows 


2 3 1 


r" Xt “1 


0i = 


^1 - l) *f ^ x 3 


to « 


dM 

X 3 X 5 5 


x <f x l ‘I) - 3 XaX 3 


the final differential equation is obtained 


a 2 

ax 2 


r~?r - 0i 2y 2 = 02 


(17) 


On assuming $g equal to a positive constant, the solution is 


02 

Ug = Cj_ sinh $ 2 . x + C 2 cosh <f)^x - — 2 

01 


( 18 ) 


The constants of integration are determined as follows: 

aug 


At x = 0, 


- 0 


therefore Cj = 0. At x = L, 


«2 " 0 
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theref ore 


Consequently 


Substituting 

yields 


Uq = 


At x = 0 , 


therefore C 


C 2 = 


0 ]_ 2 cosh 0 ]_L 


*2 "72 


02 f 30 sh 0 ^x 


01 2 \cosh 02 .L 


- 1 


dx 2 


Ug and — ■ into the third equation of equation 


d u q 0 2 /cosh 0^x \ Xi* 0 2 cosh 0jX 

^x 2 X3 0-^ 2 \cosh 0^L 


X3 cosh 02_L 


duo 0 2 sinh 0 x x X^ 0 2 sinh 0jX 0 2 x 
dx X3 0i^ cosh 02_L X3 0j_ cosh 0j_L X^ 


0 o cosh 0 t x Xj. 0 2 cosh 0 ix 02 x 2 

— — ±_ - — - + CoX + 

4 a T , a 2 ^ t ^ A 2 3 


X3 0-^ cosh 0^L X3 0^ 2 cosh 0^L 2X3 0^ c 


^0 

dx 


= 0 


= 0. At x = L, 


( 19 ) 
( 15 ) 

( 20 ) 
; ( 21 ) 
Cl* ( 22 ) 


1^ = 0 
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Therefore 


C 4 = 


0 2 


r 

Xq 02 02^ 


Xo 02.^ X3 02_^ 2X^0-^ 


Substituting into equation ( 22 ) gives 


«0 


- \?Ch ■ h) c -^ x - —,r~o - ^ - — a 2 - * 2 > (^3) 

*■301 V^i cosh ^3^1 \0q 2 2X3 0q 2 


From equations ( 19 ) and. (23), 


8u£ $2 sinh 0q x 

81 ^ cosh 0,L 

8uq ^ 02 

8x X3 ^2. 

Consequently, from equation ( 7 ), the strain is obtained as follows: 



slnh 0qx x 

cosh 0qL 01 


e x = 


0 2 

x 3 (<! 




sinh 02- x x 

cosh 02.^ 01 


(T- y -T)r 8 - i ^ < 2 *> 

dt/ 0i cosh 0jL 


and from equations (8) and (19), 


7 xy ^ 


2y 02 


d 2 ^ 


cosh 0 qx 
cosh 02_L. 


( 25 ) 
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DESCRIPTION OF TEST 


For the purpose of studying the shear-lag effect and stress-strain 
distribution on a thin plyvood sheet reinforced by longitudinal stringers 
used for the chord member of a box beam, the test specimen shown in 
figure 2 was built and tested as a simple beam loaded at midspan. 


Test Specimen 

The beam was 96 inches long, 32 inches wide, and 6 inches deep, 
open on one side and covered on the other by a plywood sheet -stringer 
combination. The box beam itself was reinforced by eight transverse wood 
stiffeners spaced 12 inches apart and four steel rods, in order to eliminate 
lateral deflection. (See fig. 2.) Five longitudinal stringers spaced 
5 inches apart were used for the plywood sheet -stringer combination. 

The side members were rectangular beams of California sugar pine, 
one with 6- by 2-inch and the other with 6- by 1-J-inch cross sections. 

Their widths were made inversely proportional to their respective moduli 
of elasticity to meet the condition of elastic symmetry for the composite 
box beam. 

The chord member was a sheet -stringer combination consisting of 
five by |-inch longitudinal stringers of California sugar pine glued 

on one face of a sheet of plywood of 0.07-inch thickness. The plywood 
sheet was composed of two 0.02-inch mahogany face plies and one 0.03-inch 
core of yellow poplar with the grain of the core and face plies at right 
angles to each other* 

The transverse stiffeners were 2- by 2-inch pieces of California 
sugar pine- The steel rods were of |-inch diameter fastened at their 

ends by screws and nuts- All component parts except the steel rods were 
glued together, with the face grain of the plywood sheet making an angle 
of k5° with the longitudinal axis of the composite box beam- 


Arrangement for Loading 

The system of loading was arranged to represent a simple beam with 
a concentrated load applied at midspan* The testing machine consists 
essentially of two pairs of jacks, individually operated, to deflect the 
ends of the beam, and of a yoke around its center connected to a lever 
balance system to measure the necessary reaction force (fig* 3)* For 
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the purpose of assuring an evenly applied load, however, one pair of 
the jacks was replaced by a stationary end support. (This alteration 
was only permissible under the assumption of very small deflection which, 
in the present case, amounts to only 0-003 percent of the length of the 
beam. ) The arrangement is shown in figures 3 and 4. 

The load was applied by operating the two screw jacks and their 
relative motion was checked by spirit levels. Steel blocks with rounded 
heads were used to approximate a line load, steel plates were used to 
avoid crushing of wood at the point of application of the load. Details 
of the loading device are shown in figure 4. 


Arrangement for Measurement of Strain 

The strains were measured by means of electrical -resistance strain 
gages connected with a portable strain indicator. Two different types 
of gage vere used: plain gages for measuring strains in only one direction, 

and rosette gages for measuring strains in three directions at 45° to 
each other. The strains were read in microinches per inch directly from 
the strain indicator* Other details of these gages are shown in refer- 
ence 5. 

The strain gages were attached to plywood and stringers along three 
different transverse sections, as shown in figures 5 and 6. This arrange- 
ment was decided upon in order to study the following characteristics: 
high shear strains at section 1 (24 in. from midspan) , appreciable shear 
and normal strains at section 2 (12 in. from midspan), and zero shear 
strain with high normal strains at section 3 (at midspan). 

Rosette gages were placed between longitudinal stringers to measure 
both normal and shear strains, whereas plain gages were located at the 
places where only normal strains were to be measured. Gages were closely 
disposed on both faces of section 2 to provide a thorough check of theory. 
However, at section 3 where zero shear strain was expected, 13 plain gages 
and only 3 rosette gages were used. Theoretically, maximum shear would 
have appeared at sections near the end of the beam, hence rosette gages 
were used at section 1 to check this phenomenon. (Local irregularities 
at sections nearer the free edge would make reasonable readings impossible-) 

The strain gages were connected to the strain indicators by means 
of a switch box to facilitate reading. (See fig. 7*) For every loading, 
corresponding strain readings of all the strain gages were taken. The 
initial zero load strain readings were checked after every series of 
loadings. The strain gages were divided into two groups and readings 
were taken simultaneously with two individual strain indicators. This 
scheme was arranged to reduce the time required for taking readings, 
thus reducing the possible influence of creep. 
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Procedure of Test 

The team was loaded by successive 200-pound increments of loads 
of: (l) from 0 to 400 pounds, (2) from 0 to 800 pounds, ( 3 ) from 0 to 

1200 pounds, and (4) from 0 to 2000 pounds. 


DISCUSSION OF FISULTS 


A very satisfactory agreement between the theoretical and experi- 
mental results for both normal and shearing strains was obtained. The 
maximum difference of strains between experimental and theoretical 
results (tables 1 and 2) was found to be about 15 percent of the experi- 
mental valuej because of the nonhomogeneity of wooden structures, this 
result is likely to be about the best obtainable. 

The comparative results were plotted as shown in figures 8 to 13 . 

By plotting the experimental normal strains for sections 2 and 3 (figs* 9 
and 10) , local shear-lag effects on the sheet between the longitudinal 
stringers were observed. This effect was to be expected, yet it was not 
taken into account in the theory. To have done so would have complicated 
the mathematics considerably. The distribution of experimental shearing 
strains resembles a cubic parabola, whereas the theoretical curve is a 
straight line (figs. 9 and 10) . Nevertheless, the agreement is very good 
so far as the maximum values are concerned. The other points on the 
curves show the theoretical results to be conservative as to magnitude 
of stresses. 

In figures 12 and 13, comparative normal strain distributions on 
stringers were plotted. In the development of the theory, it was assumed 
that the centroids of the stringers coincided with the midplane of the 
sheet. However, in practical construction this is never true. Thus a 
local bending effect of the stringers was to be expected, the result of 
which was that strains on the free face of the stringers were of a much 
lower value than those on the face to which the plywood was glued. A 
relation between them can be obtained as follows. 

If the stringer is considered as a small cantilever beam subjected 
to a horizontal load of p pounds per inch acting on one of the faces, 
as shown in figure 14, and if the vertical component of p which 
appears as a result of the bending deformation of the stringer is neglected 
by considering that the bending modulus of elasticity of the sheet is 
negligible, the normal force and the bending moment at any section can 
be given by the following expressions: 


dN x = p dx 


N x 


x 

p dx 


0 
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dM x = a p dx 

M x ~ I a p dx = a j p dx 


But. 


L xa 


-^xb 


M x a Nj 
7 + A 


Mx^ N x 


I A 

where fxa fxb are stresses. 

For constant section and constant E the ratio of strains is given 

n x 


t>y 




^b 

N x 

ab 

I 
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" A 

J pd * 

exb _ 

^xb 
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A 



Ho 

e xa 

^xa 

’ M x a 

* v 

II 

Prol 

1 

+ — 

f P 



I 

A 

I 

A 

. . 


Hence, for symnetrical sections, 


e xb 
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ab 2 
I ' A 


xa 1 

T A 


< 1 


For the particular case of rectangular stringer a - b - 
dh 3 

I = and A = dh, 


e xb 
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e xb 

Although the relation < 1 has been confirmed, the experimental 

xa e xb 1 

results shoved a large discrepancy from the relation ~2 f especially 

e xa 

for the central stringers. For design purposes, the maximum strain on 
the stringer can be assumed equal to that on the sheet- 


CONCLUSION 


The theoretical and experimental investigation of shear-lag action 
presented is based on the concept of an idealized structure in which 
the main sources of strain energy of the real structure are maintained 
while the secondary sources are neglected, and thus the idealized chord 
member has infinite rigidity in the transverse direction. 

The test results justify the adequacy of the theory for practical 
design, because the differences between theoretical and experimental 
results in the critical regions are smaller than would be expected from 
the usually irregular behavior of practical structures made of wood and 
plywood. 


Massachusetts Institute of Technology 
Cambridge, Mass., August 19^7 
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APPENDIX - EVALUATION OF PROPERTIES 
Determination of Bending Modulus of 
Elasticity of Side Members 


In order to determine the bending modulus of elasticity of the side 
members, a transverse bending test was performed (fig* 15)* The specimen 
was tested as a simple beam, supported near the ends and loaded at 
midspan. The tests were run with the specimen in normal position and 
with it turned over. Deflections at midspan were measured by means of 
a dial gage. Load-deflection curves were plotted and from them the 
bending modulus of elasticity was computed. A correction term for deflection 
due to shear was included in the formula for E* 

The average values for the modulus of elasticity were found to be. 


Beam 1 1,782*000 psi 

Beam II 2,026,000 psi 


Determination of Tension Modulus of 
Elasticity of Longitudinal Stringers 

The value of E for the longitudinal stringers was determined 
from tension tests on specimens having the dimensions shown in figure 15* 
The specimens were cut from three individual longitudinal stringers, the 
longitudinal ajces of which were parallel to wood grain. Tension loads 
were applied with the ends of the specimen attached by tightly clamped 
steel plates. Strains were measured by electrical -resistance strain 
gages with strain indicators, and Huggenberger tensometers were set on 
both sides of the test specimen for check readings. 

Load-strain curves were plotted, from the slopes cf which E was 
computed. The average value for the three specimens was found to 
be 1,1+96,000 psi. 


Determination of Modulus of 

Elasticity and Poisson’s Patio for Plywood 

Modulus of elasticity and Poisson’s ratio were determined by tension 
tests similar to that outlined for determination of E for longitudinal 
stringers. A simple bolted joint was used a z the ends of the specimen. 
Electrical-resistance strain gages were glued par-all el to the axis of 
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loading on one face and perpendicular on the opposite face of the specimen 
to measure longitudinal and lateral strains. The specimens were cut 
out of the same sheet used for the composite box beam. The sizes were 
as shown in figure 1 6 , with the longitudinal axis parallel, perpendicular, 
and at a ^5° angle with the face grains. Load-strain curves were plotted 
from which the data on modulus of elasticity and Poisson's ratio were 
computed. 

The reading from the longitudinal gage is affected by transverse 
strain and that from the lateral gage is affected by longitudinal strain. 
In order to take into account these effects, the formula for rosette 
strain gages was applied, and longitudinal and lateral strains were 
obtained from 

e 3 

6 l - e i - 55 


e i 

€ T = e 3 ‘ 55 


where 

e^ longitudinal strain-gage reading 

e^ transverse strain-gage reading 

Since is always small compared with e^> the second term could 

be neglected in the formula for and E computed directly from 

ei readings. 


p/a p/a 


However, the effect of longitudinal strain on lateral strain is of an 
appreciable amount; consequently Poisson's ratio must be computed by 




e L 



0.02 
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From the tests, 


J uv = 0.1152 
J vu = 0.1189 
E u = 875,000 

Ey = 837,000 


Determination of Modulus of Rigidity 
for Plywood 

Based on a well-known result in the theory of tending of thin 
plates (reference 6) , a method of determining modulus of rigidity was 
outlined in reference 7« The method consists in determining the deflection 
of a square plywood plate loaded at two opposite corners and supported 
at the other two. The deflections were measured at points on the diagonals 
equally distant from the center. Load-deflection curves were plotted, 
from the slopes of which the modulus of rigidity was determined from the 
formula 


_ 2 I 

2 w 

h“ 


where 

P load on each corner 

w deflection of points on diagonal relative to center of plate 

u distance from center to points on diagonal 

h thickness of plate 

The specimens used were 3" by 3"i nc k and 2^- By 2^-inch sizes and 

their edges were either parallel or perpendicular to the face grains. 
Small thin copper sheets were glued on the corners in order to eliminate 
local stress effects at load points. 


The defied tions were measured by means of dial gages- The load 
was applied by a testing machine as shown in figure 17* Effects on dial 
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reading caused by Initial curvature and gage spring force were minimized 
by obtaining two sets of readings for each specimen, one with the loads 
acting downward on two specified opposite corners, the other with the 
specimen rotated through 90° so that the loads acted upward on these 
corners. At the same time, dial gages were carefully chosen so that 
their spring forces were about equal. 

From the tests, the following results were obtained: 

3" hy 3-inch specimens 2^- by 2^-inch specimens 

G uv 140,300 psi 139,300 psi 


An average of five specimens gave a value of G uv of 139,900 psi. 


Stress -Strain Relations for Plywood 
at Any Angle to the Grain 

It is usually assumed that plywood behaves elastically as an 
orthotropic material. The axes parallel and perpendicular to the grain 
are generally taken as principal axes of strains. 

Let E u and be the moduli of elasticity 

parallel and perpendicular to the face grain} J vu 
ratios in those directions} and G uv' the modulus 

to shearing stresses parallel and perpendicular to 
respectively. 

With x and y indicating directions at 45° to the face grains, 


in the directions 
and J uv , the Poisson's 
of rigidity referred 
the axes u and v. 



J xy 


= J. 


yx 


1 _ ^ ~ ^uv _ ~ ^VU 

Guv E u Ey 

1 7 ~ J uv 1 ~ J vu 

G uv E u E v 


> 


1 _ + **uv 7 + J vu 

G xy ~ E u + E v 


J 


( 26 ) 
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A derivation of these relations can he found in references 2, 8, and 9* 

From the tests, it was found that E^y = 546,000 and 555,000, 

Jjgr = 0.520 and 0.609, and G^y = 383,000. 


COMPARATIVE VALUES BASED ON DATA FROM REFERENCE 10 
(l indicates parallel to grain; T indicates perpendicular to grainj 



Mahogany 

Yellow poplar 

E L (psi) 

1,380,000 

1,430,000 

Et/El 

0.039 

0.039 

Erji( psi) 

54,000 

52,000 

J TL 

0.552 

0. 406 

g^t/e 

0.038 

0.053 

G^ip(psi) 

53,500 

76,000 


E u = “ (t^E^ + ^ ^4 X 1380) + (3 X 523 X 1000 


= 811,000 psi 


E v = 1 (t 2 %. m + ^ I 4 X 54) + (3 x 1430)] x 1000 


= 64U,000 psi 
J vu = ~ (^ETmJLTm + ti^p^L^) 

= — - — 1(4 x 54 x 0.552) + (3 x 1430 x 0.020)] = 0.0455 
7 x 644 L 

J uv = (^ E ]^TLm + t l E Tp J LTg) 

= x 1380 x 0.0215) + (3 x 52 x 0.4067) = 0.0322 
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G = 1 
^uv t 


(^2^m + t l G p) 


= [(4 x 53-5) + (3 x 76)] x 1000 = 63,200 psi 


_1 _ 1_ 
El 


1 - J-^y 1 * J" 


Ey 4 E. 


u 


4 E„ 


vu 1 
+ 


4 G 


uv 


1 ^. 0; . 0 i g g t 1 - . 0.0^5 „ 1 4. 626 x lo~^ 

4 x 811,000 4 x 644,000 4 x 63,200 


E x = Ey = 216,000 psi 


! 1 - J- 

^uv 


UV 1 ^vu 


E-n E v 

j xv = j vx — - — 

** yX 1 1 - J U y 1 - Jy U 


x uv 


E 


u 


E„ 


_1 1 - 0.0322 _ 1 - 0. 0455 

63.2 ~ 8ll 644 

1 1 - 0.0322 + 1 - 0.0499 

63.2 + 811 + - 644 


0-712 


1 


1 + J 


uv 


1 + J 


vu 


J u 


E, : 


1 + 0.0322 
811,000 


1 + 0.0455 „ -6 

— = 2.897 x 10 ° 

844,000 


G-^ = 346,000 psi 
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Check of Strain-Energy Relation E U J UV = EyJ TU 



Vvu 

E v J vu 

Difference 

Error from 
mean value 
(percent) 

Experimental 

value 

100.6 x 10 3 

99-6 x 10 3 

1.0 x 10 3 

1.0 

Reference 10 

26.1 

29.3 

3.2 

11.6 


The average or standard data given in reference 10 do not conform 
to the requirements of the strain-energy relation as veil as the data 
from a limited number of tests on the material used in this beam. An 
effect of this sort would be expected, perhaps, but the actual discrepancy 
when the data of reference 10 are used is greater than would be anticipated. 
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SUMMARY OF ELASTIC PROPERTIES OF WOOD MEMBERS 



From 

experiment 

Computed 

from 

experiment 

Computed 
from 
data of 
reference 10 

| Difference from 
data of refer- 
ence 10 
( percent) 

Side members (California sugar pine) 

E^psi) 

1782 x in 3 


1144 x lo3 

56.0 

E 2 (psi) 

2026 X 10 3 


1144 x 10 3 

72.0 

Stringers (California sugar pine) 

Egt(psi) 

1496 X 10 3 


3 

1144 x 10 

30.6 


Plywood ( mahogany j yellow poplar) 


E u (p8i) 

875 x 10 3 


811 X 10 3 

8-3 

Ey ( psi) 

837 x 10 3 


644 x 10 3 

14.4 

*^uv 

0.115 


O.0322 




J vu 

0.119 


0.0455 




G UY (p3i) 

139-5 x 10 3 







Exy(psi) 

550.5 x 10 3 

433 x 10 3 

216 x 10 3 

155 

J xy 

0.565 

0.553 

0.712 

-6.6 

G xjr (psi) 


383 x 10 3 

346 x 10 3 

10-7 



Because of the length of the tables recording the observed strains 
and of the computations of the theoretical strains, test data and 
computations are not recorded here. The reduction of the test readings 
was done by standard methods, by using the properties of wood obtained 
from experiment rather than those from reference 10. The theoretical 
values were computed from the procedure presented in the first part of 
this report, again by using the properties of the wood in this particular 
beam. 
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The curves of figures 8 to 13 show the agreement between theoretical 
and experimental strains when the values resulting from the computations 
are plotted. The agreement is, in general, satisfactory evidence of the 
accuracy of the computations which are not included here. 
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TABLE 1 


<_o 

4 =" 


THEORETICAL STRAINS 

[ior 1000 It at midspan} for 500 lb at tip] 


X 

(in 0 from 
free edge) 


Sheet 

Stringers 

( in .) 

€ x 

^xy 

€ st 

(top) 

e st 

(hottom) 

12 

0 

5 

10 

12.1875 

14 

64.01 x lo -6 

680 85 

83.36 

92.74 

101.93 

-552.68 x 10-8 
-634.87 

64.01 x 10" 6 

68085 

83-36 

-6 

32.O x 10 

34.42 

41.68 





22 

(Section l) 

0 

5 

10 

12.1875 

14 

110.27 

120.33 

150.49 

170.01 

189.10 

- 485.43 

-557.62 

110.27 

120.33 

150.49 

55.13 

60.16 

75.25 





34 

(Section 2) 

0 

5 

10 

12.1875 

14 

146.19 

165.80 

224.60 

262.66 

299.88 

- 316.37 

-363.42 

146.19 

165.80 

224.60 

73.09 

82.90 

112.30 





46 

(Section 3) 

0 

5 

10 

12.1875 

14 

141.920 

177.815 

285.50 

355-18 

423.34 

0 

0 

141.920 

177.82 

285.50 

70.96 

88.91 

142.75 
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TABLE 2 

EXPERIMENTAL STRAINS 


|?or 1000 11 at midspanj for 500 11 at tipj 


Section 

y 

(in.) 

Sleet 

Stringers 

€ 

X 

6 

y 

7 

xy 

6 at 
(top) 

e st 

(bottoni) 

1 

-12.1875 

- 7.5 

- 2.5 

182 . l 4 X 10 -6 

1 V 7.35 

145.31 

-11.80 X 10 ~ 6 

-25.26 

-22.23 

- 496.40 x 10" 6 
-201.03 
-15.28 









2 

-15 

-12.1875 

-10.0 

- 7-5 

-5 

-2.5 

0 

2.5 
5.0 

7.5 

10.0 

13.00 
15.125 

354 - 

229.64 

270 

190.35 

195 

138.76 

164 

137.26 

180 

191.30 

264 

240.49 

353 

-70.36 

- 294.42 

354 X 10 -6 





270 

176 x lO- 6 

-25.67 

-H7.53 

195 


10.16 

- 19.97 


164 

-26 

10.19 

8.05 

180 

50 

-34.69 

117.50 

264 

81 

-35.69 

302.40 

353 





3 

-15 

-12.1875 

-10 

- 7.5 
-5 

-2.5 

0 

5 

10 

15.125 

479 

304.53 

300 

205.15 

202 

152.23 

160 

193 

219 

472 



479 


-73.85 

- 16.33 


300 

183 

-22.08 

13.13 

202 

54 

32.62 

4.02 

160 

193 

291 

472 

-17 










00 

ui 
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Figure 1.- Symbols used in the theory. 
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Figure 2.- Test beam. 
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Figure 3.- Beam ready for test. 



Figure 4.- Detail of loading device 





Section 1 
(top) 

24 inches from 
midspan 


Section 2 
(top) 

12 inches from 
midspan 


Section 2 
(bottom) 


Section 3 
(top) 

at midspan 


Section 3 
(bottom) 


XST Tkzi YVW 



12 H 3 H 4 

1 I |t 



% *■ “ I ir - | 




/ 6 ~ 






Figure 5.- Location of strain gages. 
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(a) Top. 



(b) Bottom. 


Figure 6.- Location of strain gages 
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Figure 7.- Strain indicators and switch boxes. 
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Figure 8.- Comparison of theoretical and experimental strains. Section 1. 

x = 22 inches. 


in. /in. 
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x = 34 inches. 
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x = 46 inches. 


, in. /in. 


cn 

O 



Tip 10 20 30 40 Root (or midspan) 

x, in. 

Figure 11.- Comparison of theoretical and experimental spanwise strain distributions. 
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Top 

~W 

Bottom 


m 



Figure 12.- Comparison of theoretical and experimental strains in 
stringers. Section 2. x = 34 inches. 
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Top 

s ^ m 

Bottom 



Figure 13.- Comparison of theoretical and experimental strains in 
stringers. Section 3. x = 46 inches. 





Figure 14.- Notation for simplified plate -stringer. 
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(a) Test specimens for determining G. 



(b) Detail of loading device for determination of G. 


Figure 17.- Test setup and specimens for determining shear modulus G. 
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